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CN ' Abstract 

Nielsen [3] recently asked the following question: "What is the minimal size quantum 
circuit required to exactly implement a specified n-qubit unitary operation U, without 
the use of ancilla qubits?" Nielsen was able to prove that a lower bound on the minimal 
ON ■ size circuit is provided by the length of the geodesic between the identity I and U, where 

the length is defined by a suitable Finsler metric on SU(2 n ). We prove that the minimum 
circuit size that simulates U is in linear relation with the geodesic length and simulation 
| parameters, for the given Finsler structure F. As a corollary we prove the highest lower 

bound of 0{^d 2 Fp {I, U)L Fp (I, U)) and the lowest upper bound of n(n 4 df, (J, 17)), for the 
standard simulation technique. Therefore, our results show that by standard simulation 
one can not expect a better then n 2 times improvement in the upper bound over the 
result from Nielsen, Dowling, Gu and Doherty [4]. Moreover, our equivalence result can 
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■ be applied to the arbitrary path on the manifold including the one that is generated 

*^ adiabatically. 

1 Introduction 

& : 

Quantum computation is inherently a process of continuous evolution of quantum states that 
has the potential to fundamentally change the notion of feasibly tractable computation. Only 
recently did researchers start to think how notions from the differential geometry [13] can be 
^ ■ used to represent this process. Instead, quantum circuits, as an inherently discrete notion 

of computation, have been most commonly used to represent this continuous process. Any 
quantum operation, a unitary matrix U G SU(2 n ), is an element of a Lie group, and a point on 
the IA = SU(2 n ) manifold, whose tangent bundle can be endowed with the Finsler structure 
F, that effectively provides a measure of length for any path on the Finsler manifold (U, F). 
In particular, the paths that we are interested in are geodesies. These are locally and, under 
certain conditions, globally minimal length paths between any two points on the manifold. 
They are of particular interest, because if closely simulated they can give the smallest circuit 
complexity for the given unitary U. 

The aim in the approach that we take here is to tackle the question about the complexity 
of the circuit necessary for the simulation of an arbitrary unitary gate. As the length of the 
geodesic for the particular unitary is its intrinsic property, ideally one would succeed in finding 
the minimum number of circuits necessary to implement the unitary by simulating exactly 
its geodesic. Therefore, the hope here is to learn about the circuit lower bounds by basically 
transforming the hard combinatorial optimization problems over large sets to the problems in 
continuous domain that can be solved with tools of differential geometry and the calculus of 
variations. 
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One of the first results that had the flavor of this transform was introduced by Mochon 
[5], who proved that in the discrete model and the analogous continuous model, i.e. the 
Hamiltonian oracle model, Oracle interrogation, the problem of computing XOR and Grover 
search have the same complexity. Moreover, Nielsen [3] and subsequently Nielsen, Dowling, Gu 
and Doherty [4] proved that for particularly chosen metric there is a polynomial equivalence 
between the geodesic length and number of gates necessary for the simulation. The lower 
bound for the minimum number of gates necessary for the simulation has been determined for 
exact simulation and the upper bound has been determined for the arbitrary precision. The 
metric chosen in ]4J penalizes all those directions on the manifold that are not easily simulated 
by local gates, so that coefficients for stabilizer elements of Hamming weight greater than two 
bear high cost, i.e. have longer paths. 

In this paper we prove the stronger result and show the exact upper and lower bound 
that determine the equivalence between the minimal number of gates in the standard circuit 
simulation and the length of the geodesic. Both upper and lower bound are determined by 
the simulation parameters and, of course, the length of the geodesic. 

We consider the two cases. First: the simulation of the geodesic with set of gates Q that 
is exactly universal, and the case with approximately universal set of gates. With the exactly 
universal set of gates for any point xq 6W there is a gate in a set Q by which we can simulate 
exactly any point x\ E IA in the ball of radius r centered around xq; we denote this ball 
as B^ Q (r). Under that assumption, we prove that the number of gates in the simulation of 
a geodesic may be upper and lower bounded by a linear factor in the length of geodesic and 
simulation parameters. 

When the set of gates, Q e , is approximately universal, a single gate from this set can 
simulate the points in B^ (r) only with some finite precision e, and that will necessarily mean 
that the circuit that simulates the geodesic is doing so along the path that is not shorter than 
the actual geodesic, for that very point. 

Our aim here is not to elaborate on the algorithm for the geodesic simulation but rather 
to prove the bounds that optimal simulation can achieve. We say optimal, because the set 
of gates Q that we first consider is much more powerful than any local and universal set of 
gates. Therefore the result that we present is the optimal result about complexity equivalence 
between discrete and continuous notions of computation. In particular, for the standard 
simulation model described in [4], we derive the highest lower bound and the lowest upper 
bound in the minimal circuit complexity that one can hope to achieve with the simulation of 
a geodesic. 

2 Preliminaries 

A quantum operation U G U is a point on the manifold hi = SU(2 n ) at some distance 
from identity / G U. The distance considered is the integral distance that is determined 
by the structure used on the manifold. In general that structure may be more general than 
Riemannian, i.e. it is called the Finsler structure F(x, y). The restriction of a Finsler structure 
F to any specific tangent space T x hi with the origin at the point x G U is called Minkowski 
norm on TJJ. The second argument of the structure F(x,y) is the velocity and its definition 
follows. Therefore a Finsler structure is basically a family of the smoothly varying Minkowski 
norms, one for each tangent space. 

The defining properties of a non-negative real-valued structure F(x,y) on M 4 " -1 are as 
follows: 
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(1) it is C°° anywhere on R 4 " -1 except at y = 0; 

(2) it is positive homogeneous, i.e. F(x, Ay) = XF(x, y) for A > 0; 

(3) the (4 n — l)x(4 n — 1) matrix -§^--^p[\F 2 ] is positive definite unless y = 0. As a conse- 
quence, one can derive positivity and triangular equality of Minkowski norms [6]. The 
structure F(x,y) is usually denoted simply as F(y). 

For any a, b G M+ we say that a map a : [a, 6] — > U is a piecewise C°° curve with velocity 
y = II = Si !Iir~tE~ e T a (t)U- The integral length of the curve er, L(cr), is defined as: 

L F (a) = J F(a,^)dt. (1) 

Since we are usually interested in minimum length curves for xq, x\ € U, we denote by T(xq, Xi) 
the collection of all piecewise C°° curves a : [a, 6] — > U such that cr(a) = xq and cr(6) = xi. 
Similarly, the integral distance is defined as a map d F ■ U x U — > [0, oo): 

dp(xo, xi) = inf Lp(cr) (2) 

r(a; ,a;i) 

Using these definitions, one can show that the Finsler manifold (W, c^f) satisfies the two axioms 
of a metric space: (1) positivity: dp(xo, x\) > 0, where equality holds if and only if xq = x\ and 
(2) the triangular inequality: dp{xo,X2) < dp(xo,xi) + dp{x\, x^)- In general, the symmetric 
property of a distance does not need to hold, and therefore dp(xo,xi) ^ dp(xi,xo). 

3 Distortion Lemma 

To establish the equivalence result, we introduce in this section the main tool of our analysis. 
The intuitive idea on which we build our results relies on the relation between the distances 
on the manifold and the distances on the tangent space of the manifold. While the former are 
introduced by the unitary gates and their complexity, the latter are defined by the appropri- 
ately defined distances between the Hamiltonians of gates used in the simulation. This will 
be proven useful in the sections below. 

The lemma that follows is a slightly stronger result of a well-known and very useful fact 
from the differential geometry. Again, it relates distances on the manifold with the minimum 
and maximum distortion of the Euclidian norm on the tangent space over the compact set. 
Interested reader are encouraged to consult [6j, an excellent and very elaborate reference on 
this subject. 

Lemma 1 (Distortion Lemma). Let (14, F) be a Finsler manifold, and for any point i£W 
let tp : P x — > M 4 ™ -1 be the local coordinate system diffeomorphism of a compact set P x onto 
an open ball o/M 4 " -1 , such that ip(x) = 0. Then for a given xq, x\ G P x , and any Finsler 
metric F(x,y), there exist a constant minimum m > and a constant maximum 9Jt > 1 such 
that the following relation is true: 

m\tp(xi) - <p(xo)\ < L F (x ,x 1 ) < 3Jt\tp(xi) - if(x )\ . (3) 

Here \(p(x\) — <p(xo)\ denotes the Euclidean length of the 4 n — 1 dimensional vector in the 
tangent space. 
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Proof: We first note that a compact set P x for which tp(x) = always exists. This is 
true because, given a local coordinate system tp : Q — > 1R 4 " -1 and x G Q for which tp(x) = 0, 
we can choose P x to be a closure of the preimage of p^ 1 (B 4 _1 (r)) for some r > 0. By 

B^ n ~ l (r) = {v G 1R 4 " -1 : |«| = yJ2i v 1 < r l we denote the ball of radius r in the tangent 
space whose closure is a subset of tp(Q). 

Next we note that, for tangent vector y = Y^iVi'&r = W ^ T X U, the ratio between 

Minkowski norm F(x,y) and x-dependent Euclidean norm \y\ := \J^2i yf for the basis {^-} 
is well defined for y ^ 0. Since both norms are positive continuous functions over the compact 

sets their quotient is also a positive continuous function. Therefore the quotient's minimum 

F(y) 

>L anu cll^ v\j\jll ^uaiuivc. \J - — III _^ 

y G T X U and all x £ P x : 



m and maximum 971 exist and are both positive: < m < -r*p- < 971. In other words, for all 



m\y\<F(x,y)<m\y\. (4) 

Now we can prove the right hand side (RHS) of inequality (j3|) by choosing the path a G P 
that maps under (p to a line segment. In that case we can write: 

/•ti rti 
L F (x ,x 1 )= F(a')dt<m \a'\dt = Tl\tp(xi)-<p(x )\ , (5) 
J to J to 

where a' = denotes the velocity field of a path a. 

To prove the left hand side of inequality we first show that a must be contained in P x . 
The proof is by contradiction as follows. 

Choose ro < m ^, gn r and eo = 9Jtro and Pq = Lp~ 1 [B n {ro)] C P x . Let a : [to,ti] — » W 
be a piecewise C°° curve such that a (to) = xo and a(t±) = x\ for xq, x\ G Po- If Lp(a) < 
dp(xo,xi) + eo then the curve a is certainly contained in P x , and since by equation © 
dp(xo,xi) < 29Jtro we have by assumption that Lp(a) < 39Jtro. Now if we suppose that a is 
not contained in P x , and let to < t* < t% be the first instance where a reaches the boundary 
dP x , at the point q = cr(t*), so that \ip(q)\ = r, then: 



L F (a) > L F (a M ) = I F(a')dt > m f \a'\dt 

J tg J to 



>m\(p(q) -<p(x )\ >m(r-r ) . (6) 



But the length of this curve would in fact be longer then the maximum possible length of 
39Jtro < m(r — ro), since by assumption we are assured that for m > and 9JT > 1 it is true 
that ro < m4 ^ OT r. Therefore a must be contained in P x . 

The proof of the left hand side of inequality ([3|) follows by the same arguments as were 
used to prove (JBJjD 

Given the distortion lemma for the length of the path for any two points that belong to 
the compact set, we can easily derive a simillar result that is valid for the shortest distances. 

Corollary 1. For a Finsler manifold (U,F), and any point x G IA, let tp : P x — > IR 4 " -1 be the 
local coordinate system diffeomorphism of a compact set P x onto an open ball o/lR 4 "" 1 , such 
that tp(x) = 0. Then, for a given xq, x\ G P x and any Finsler metric F(x,y) there exist a 
constant minimum m > and a constant maximum 9JT > 1 such that the following relation is 
true: 



m\tp(xi) - <p(x )\ < d F (x ,xi) < Tl\tp(x t ) - tp(x )\ . (7) 
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Proof: We only need to verify the left hand side of inequality ([3]) is still true for minimal 
length curves. By definition of metric distance, for < e < eo, two points xq,x\ € Po can be 
joined by a piecewise C°° curve a : [io>*i] U with integral length: 

L F (a) < d F (x ,xi) + e . (8) 

By previous arguments, a must lie in P x , end by similar calculations we find that: 

m\(p(xi) - (p(x )\ < L F {a) < d F (x ,xx) + e , (9) 

Letting e — > proves the desired result. □ 

Lemma ([1]) and Corollary ([1]) allow us to bound the lengths on the manifold to the Euclid- 
ian lengths on the tangent space. For Euclidean coordinates in our tangent space we will have 
the coefficients in the decomposition of the gate Hamiltonian matrix in terms of the generalized 
Pauli matrices, n times tensored two dimensional matrices from the set {I,X,Y,Z}. 

For example, in the context of simulation, points xq, x\ G P x and an open set P x are 
chosen such that they correspond to the end points in the simulation by a single gate. More- 
over, we can construct a local coordinate system on the Lie group SU(2 n ) which is a Lie 
algebra su(2 n ). For the origin x s € U, define a pull back map (p^ 1 : IR 4 ™" 1 — > U, so that 
x s+ i = exp~ lVs+1 ' (T x s = exp~ ltp ( Xs+1 ^' (T x s , where a denotes the coordinate basis, i.e. (4 n — 1)- 
component vector whose entries are the generalized Pauli matrices. 

For the particularly chosen metric, as in [4J, F p (x ,y) = F p (y) = ^Ya=i vf + P 2 X^i vh 

where k = 9 ^ n 2 n ^ +3n, which introduces a penalty p for the subset of Hamiltonian coordinates 
in the tangent space, so we have: 

\y\ < F p (y) < p\y\. 

Since this relation is true on any compact set, by Corollary (pQ) we have: 

\(p(x s+ i) - <p(x 8 )\ < d Fp (x Sl x s+ i) < p\ip(x s+1 ) - <p(x s )\ . (10) 

It is important to note that in our analysis constants m and 9Jt do not depend on the compact 
set within which each gate is applied, and they are basically the property of the metric. This 
property might not be true in general for some other Finsler structures, but for our purposes 
here this assumption is very plausible. 

4 Equivalence Result 

For the sake of consistency and easier understanding, we follow the notation from [3] and 
denote with mg the minimum number of gates, for a given set Q, needed to implement an 
arbitrary unitary U ZLlA. Moreover, in this section we assume that the geodesic is simulated by 
sequential application of the gates from the set Q, and that by using a single gate from the set Q 
we can simulate exactly any other point in the ball B x (e), i.e. which is the e— neighborhood 
around the initial condition at the point xq. This is an unrealistic scenario, since the set 
of gates would need to be infinite and non-local. Hence, we relax it in the next section. 
However, for the purpose of exact simulation of the geodesic it is an important tool. Clearly 
consideration of the set of gates defined in this way, as we shall see, is the best we can possibly 
hope for, and thus the bounds achieved by the set Q are optimal, i.e. they determine the 
bounds achievable by any other set of gates that is less powerful. 
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For any gate used in a simulation we assign a gate index, so that eventually the index set 
is s = {0, 1, 2, mg — 1} for every gate in the simulation. Moreover, by a(t) : [0,mg] — > U 
we denote a minimal geodesic between cx(0) = /, a(mg) = U. 

Note that, since U is a compact manifold with Finsler structure, all forward and backward 
Cauchy sequences with respect to d must converge on U. More precisely, compact Finsler 
spaces are automatically both forward complete and backward complete. This fact holds 
regardless of whether the Finsler structure is absolutely homogeneous or only positively ho- 
mogeneous. Therefore, any two points on the manifold can be connected by a minimizing 
geodesic, as that property itself is a sufficient condition for the Hopf-Rinow theorem [UJ. 

Theorem 1. Let d F (I,U) denote a length of a geodesic between I and U £ SU(2 n ). For any 
simulation index set s = {0, 1,2, ...,mg — 1} let P Xs Gil be an open set on the manifold that 
contains a segment of minimizing geodesic a s (t) : [s,s + l] — > that is simulated exactly by a 
single gate. Moreover, let P Xs be mapped by tp diffeomorphically onto an open ball in IR 4 ™ -1 , 
so that p s = |</?(x s +i) — (f(x s )\ is the Euclidean length of the image of the geodesic segment 
a s . If we denote p sup = sup s p s and p ini = inf s p s , then the following relation holds: 

WJl <mg< tiMl, 

Proof: For any segment gate index from set s, by the Corollary (|TJ) we see that: 

mp s < d F (x s ,x s+1 ) < Wlp s , 

Summing over all segments of minimizing geodesic Y2T=o 1 d-F^Xs, %s+i) = dp(I,U), and 
taking into account that YlT=o 1 Ps — m gP an d m g^f < YlT^o 1 /T> ^ ^ s eas y t° see that: 

m g -l 

mgmp in{ < d F (I,U) = ^2 d F (x s ,x s+1 ) < mgMp sup , (12) 

s=0 

which gives the desired result by rearranging the variables. □ 

Note that the above theorem is derived in terms of bounds of the Euclidean distances 
in the tangent space. One may take a different path though, as for example Nielsen in [3], 
by deriving the result for the lower bound in terms of the lengths of the geodesic segments 
simulated by the single gate: d F (x s , x s+ i) < /3 sup . From the following theorem, one can 
reproduce the result derived by Nielsen as a special case when f3 snp = 1. 

Theorem 2. Let d F (I, U) denote a length of a geodesic between I and U S SU(2 n ). For any 
simulation index set s = {0, 1,2, ...,mg — 1} let P Xs be an open set on the manifold that 
contains a segment of minimizing geodesic cr s (t) : [s, s + 1] — > U, that is simulated exactly 
by a single gate. Moreover, let P Xs be mapped by <p diffeomorphically onto an open ball in 
R 4 ™ -1 , so that p s = \lp(x s+ \) — ip(x s )\ is an image of the bounded length geodesic segment 
d F {x s , x s+ i) < (3 S . If we denote f3 sup = sup s (3 S and /3 inf = inf s /3 S , then the following relation 
holds: 

^Ml < ms < «MLM . {13) 

Psup m Pinf 

Proof: Following along the lines of Theorem (P): 

R 

m m = mPs ~ ^( Xs ' Xs +i) - Ps = Wp s , 
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Summing over all segments of minimizing geodesic J2s=o dp(x s , x s+ i) = dp(I,U), and 
taking into account that rng(3 ini < YlT=o @s — m G@sup- 

^rng(3 ini < d F (I, U) < mgf3 sup , (14) 

which gives the stated result. □ 

Equations (fTTI) and (fTil) establish the tightest possible equivalence between the minimal 
number of gates in the circuit and geodesic length as a function of the simulation parameters. 
Again, the simulation parameters may be defined in terms of distances traversed with the 
single gate on the manifold or in terms of the Euclidean distances between the initial and 
final coefficients in the generalized Pauli expansion of the gate Hamiltonian. Even though the 
above results give no indication as to how to implement the simulation, they do provide us 
the best bounds we currently have and give us an estimate to the quality of the simulation 
provided that one knows the simulation parameters. However, the above results can be applied 
to the arbitrary paths on the manifold including those that are generated adiabatically. In 
particular, it would be very interesting to compare the results for bounds of circuit size 
obtained by geometric techniques with the equivalence results obtained in pQ. 



5 Approximate simulation 

In this section we reformulate the bounds for the standard circuit simulation procedure where 
the set of gates used consists solely of the single and two qubit gates, which are applied 
sequentially. Since the exact simulation of arbitrary unitary gate by single and two qubit gates 
demands an exponential number of gates, almost all unitaries simulated by the polynomial 
number of gates will be simulated approximately. 

In particular, we consider two paths. Let the first be dp p (I,U), denoting the length of 
the geodesic simulated exactly with the set of gates from Q with respect to the Finsler metric 
F p , and let the second one Lp p (I, U) be the minimum length path for the exact simulation 
of U by the set of gates from Q2. Here we denote by G2 the set of unitary gates whose time 
independent Hamiltonians have Hamming weight not greater than two. 

Note that the length Lp p (I,U) has nothing to do with dp p (I,U), as Lp p (I,U) is com- 
pletely determined by the simulation, and almost everywhere does not simulate the geodesic 
dF p (I,U). 

Corollary 2. Let U be the approximation of the unitary operation U that is simulated by the 
one and two qubit gates. Then the lower bound on the minimum circuit size fhg 2 is at most 
0(—dp p (I, U)Lp p (I, U)), and the upper bound on mg 2 is at least £~l(n 4 d F ^(I , U)) . 

Proof: The three step standard simulation of arbitrary U = e -4 ^*)* is elaborated in 
detail by Nielsen, Dowling, Gu and Doherty in [3J. The procedure can be sketched as follows: 

(1) the time variable Hamiltonian H(t) is substituted by projected the Hamiltonian Hp{t) 
that is formed by deleting all o~i for i > k, i.e. all three- and more-body terms in the 
Pauli expansion of H(t) = Ya=i + Sjyi Vi a ^ where k = 9(n 2 ~ w ^ + 3n; 

(2) the evolution due to H p (t) is broken up into many small intervals, each of length A, 
over which the time-dependent Hamiltonian H p (t) is accurately simulated by a constant 
mean Hamiltonian Hp, 
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(3) the mean Hamiltonian that has k terms in the Pauli expansion with coefficients 
\yi\ < 1 is simulated with a standard simulation technique [T2] using one and two qubit 
gates. 

The reader is encouraged to see [3] for full detail of the approximation result. 

For the above procedure, since SU(2 n ) is compact and simply connected, there exists 
a path L Fp (I,U) that is exactly synthesized with the gates in the simulation. By exactly 
simulated we mean that the end points of each gate in the simulation lie precisely on the path 
of length L Fp (I, U). Clearly, the length of L Fp (I, U) > d Fp (I, U). 

Now we bound length of the path segments, Lp p (i s ,x s+ i), for each of rhg 2 gates in the 
simulation. Since there exists a compact set P Xs , such that end points x s ,x s+ i £ P Xs , that 
maps diffeomorphically to the local coordinate system, we can use Lemma (TTJ and its corol- 
laries. Corroborating the arguments used to derive equation (I10p . over the compact set 
P Xg , the Finsler structure, i.e. the Minkowski norm for the Pauli expansion of H(t), is 

F p (x s ,y s ) = yX^Ll yf + P 2 vf- Its minimum and maximum distortion over the com- 
pact set P Xg are: \y s \ < F p (y s ) < p\y s \- Therefore, by the Lemma ([1]) 

\ip(x s+ i) - (p(x s )\ < L Fp (x s ,x s+1 ) < p\(p(x s+ i) - (f(x s )\. 

The same is true for any other segment in the simulation, and hence: 

TOg-l 

i~n>S 2 PM < L Fp (I,U) = l f p (xs,x s+ i) < mg 2 pp^ up (15) 

where = inf, \ip(x s+ x) - <p(x s )\, and p^ p = sup s \<p(x s+ x) - (f(x s )\. Note that we can 
always choose the s-th gate local coordinate system so that f(x s ) = 0. 

Finally, in the three-step simulation summarized above, gates at the third stage simulate 
the time invariant Hamiltonian Hp for the segment A, with coordinates \yi\ < 1. More 
precisely, the s-th gate simulates the neighborhood around x s : x s +i = e~ tLp ^ Xs+1 ^'' 7 x s = 
e -«?/ s cr s A Here a G Q denotes stabilizer basis on n qubits, a s € G2, and A 2 is the simulation 
time for every gate. If we choose A = &((n 2 d Fp (I, t/)) -1 ), as in [3], then for \y s \ < 1 we see 
that p s = \y s A 2 \ = @((n 4: d' Fp (I, U))^ 1 ). Finally, using equation (fT5|) we establish that the 

simulation with rhg 2 gates has the upper bound ®{n A d 2 F (I,U)L Fp (I,U)) > Q(n 4 d F (I,U)). 

By similar arguments, for the lower bound, we get 0(^ p -d Fp (I, U)L Fp (I, U)). □ 

6 Conclusion 

The Distortion Lemma and its corollary provide a general tool for relating distances on the 
manifold with distances on the tangent space. In this paper we have derived a generalized 
linear bounds for the exact simulation of any path on the manifold, in terms of the minimum 
circuit size and the simulation parameters. 

The equivalence between the path on the manifold and circuit size still persists in the case 
of approximate simulation, provided that the simulation parameters have the appropriate 
scaling. However, one can not expect better than n 2 times improvement in the minimum 
circuit size upper bound over the result for standard circuit simulation derived by Nielsen, 
Dowling, Gu and Doherty [3]. 

Moreover, if one defines a metric on the manifold that penalizes the hard-to-simulate 
directions on the tangent space with high cost, that cost is, in effect, translated to the increased 
ratio between upper and lower bound in minimum circuit size. 
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